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Abstract 

The interaction of the partonic fluctuation of the virtual photon in deep inelastic scat- 
tering with soft color fields describing the hadron is treated in an eikonal approximation. 
It is known that, in this approach, the small-x limit of the leading-order gluon distri- 
bution xg{x, Q^) is a constant characterizing the averaged local field strength in the 
target. Matching the next-to-leading order calculation in this semiclassical framework 
with the one-loop parton model result, we obtain the next-to-leading order contribution 
to xg{x, Q"^). It shows a ln(l/x) enhancement at small x and is sensitive to the large dis- 
tance structure of the target. The final expression is a simple integral over non-Abelian 
eikonal factors measuring the target color field. We derive a quantitative relation between 
the short- distance cutoff of this integral and the scale of the gluon distribution function 
in the MS scheme. Our calculation demonstrates that higher order contributions can be 
systematically included in the semiclassical approach. 



1 Introduction 



The small-x limit of deep inelastic scattering (DIS) remains one of the most challenging 
problems on the interface of perturbative and non-perturbative QCD. It can be expected 
that both the perturbative physics described by QCD resummation techniques [0,0 and 
the non-perturbative soft dynamics underlying the growth of total hadronic cross sec- 
tions are important for a complete picture of the small-x limit of structure functions. 
In this paper, we systematically relate models for the soft color fields, which we consider 
to be a promising tool for the characterization of the hadron, to the hadronic gluon 
distribution, which is the basic object of the perturbative treatment. 

We are interested in a region of x and where, although x is very small, the 
concept of parton distributions and the conventional DGLAP evolution equations |l| 
are applicable. For the description of the DIS process, the gluon distribution g{x,ii^), 
which dominates in the small- a; region, has to be given at some scale fi^ < Q^. Thus, the 
problem of the non-perturbative structure of the target hadron has to be addressed. 

Our basic starting point is the idea of describing high-energy processes in QCD by 
studying the eikonalized interaction of energetic partons with soft color fields 0,^. In 
this approach, which has been used for the description of DIS in Q, soft color fields in 
the center-of-mass frame of the collision mediate the interaction between the energetic 
projectile and target quarks. Progress towards the description of the energy dependence 
of cross sections in the framework of soft color field dynamics was reported, e.g., in 0. 

A closely related approach to both inclusive and diffractive DIS, which treats the 
target proton as a superposition of soft color fields, was developed in P,^. In fact, this 
semiclassical approach reproduces the treatment of if the dynamics of the underlying 
color field is modelled on the basis of the stochastic vacuum |jlO[ and a phenomenolog- 
ical ansatz for the constituent- quark wave function of the proton is made. For different 



approaches to the description of the target color fields see, e.g., [|TT|. However, in the 
following we are completely general and use no model-specific features of the target color 
field configurations or, equivalently, the wave functional of the proton. 

In the present paper, the gluon distribution g{x, /i^) at x -C 1 and /i^ ^ (where A 
is a soft hadronic scale) is calculated for a target given by soft color field configurations. 



Here 'soft' means that all momentum components of the field are 0{A). Following |]T2 
a scalar 'photon', coupled directly to the gluon field, is used as a convenient theoretical 
tool for extracting the gluon distribution. The leading order calculation gives a constant 
for X g{x, /i^), which is a measure for the averaged gluon field strength in the target. This 



is in agreement with the seminal paper of Mueller |jT2[ and with ||^ , the spirit of which 
we follow closely. 

In our opinion, it is crucial, both from a theoretical and a phenomenological perspec- 
tive, to demonstrate the viability of the approach at higher orders. However, already at 
next-to-leading order the gluon distribution is scheme dependent and a careful matching 
of the partonic calculation (we use the MS scheme) and the semiclassical calculation is 
necessary to obtain an unambiguous result. With this result, we treat problems that were 



not addressed in the one-loop calculations of [jT2|,|T^, where the scheme dependence was 
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not discussed. In fact, the problems of regularization and scheme dependence arise imme- 



diately if one attempts to translate the one-loop, unintegrated gluon distribution of |13 



into a correction to the leading term in the spirit of (see [|14| for a comparison of the 



results of and ||13[ in the case of the quark distribution). As emphasized in where 
diffractive and inclusive quark and gluon distributions were calculated in the semiclas- 
sical approach, the inclusive gluon distribution dominates the small-x region. Therefore, 
we expect that our next-to-leading order result is the dominant correction relevant as 
input for the next-to-leading order DGLAP evolution. 

In our approach, the most intricate part is the semiclassical calculation at next-to- 
leading order. Working in Feynman gauge, we employ the optical theorem and calculate 
the forward scattering amplitude. In the high-energy limit, certain diagrams can be 
dropped. The remaining contribution is given in the form of a two-gluon production 
cross section. In this way, the identification with the parton model result becomes simple 
since the dangerous high-mass region, where the semiclassical approximation fails, cancels 
explicitly. 

Let us note that, for a soft hadron governed by the single scale A, the perturbative 
expansion of the gluon distribution makes no sense. However, following ideas of (see 



also the recent calculations of [0), we can always assume that we are dealing with a 
very large target, in which case the gluon distribution becomes calculable without losing 
the interest of being genuinely non-perturbative in its origin (see Jl^ for a discussion of 



the new hard scale in a framework close to the present paper). It remains to be seen in 
how far this large-target approach will allow for a description of the qualitative features 
in the realistic proton case. 

The paper is organized as follows. In Sect. § the scattering of a scalar photon off 
the target color field is calculated in the parton model and the semiclassical approach. 
The comparison of the two results gives rise to the leading-order semiclassical expres- 
sion for the gluon distribution. In Sects. |^ and ^ the semiclassical and parton model 
calculations, respectively, are carried out at next-to-leading order. The extraction of the 
next-to-leading order contribution to the gluon distribution from the comparison of the 
semiclassical and the parton model results is the subject of Sect. ||. Section || contains 
our conclusions, and a number of technical details of the calculations are outlined in 
Appendices 0-0. 



2 The leading-order result 

In the following analysis we use a scalar 'photon' coupled directly to the gluon field as 
a convenient theoretical tool for extracting the gluon distribution \T2\. To be precise, 
the real 'photon' field x couples to the field strength tensor Ff^^^ via the interaction 
Lagrangian 

Ci = ~xtTF,,F^\ (1) 

The leading-order amplitude for the scattering of the 'photon' off a classical color 
field is given by the diagram in Fig. |l|. 
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Figure [l|: The process x — > in an external color field with one gluon exchange. 

Let q and k be the momenta of the incoming virtual 'photon' (g^ = — Q^) and 
the outgoing gluon respectively. We define the light-cone components of a vector p by 
P± = Po ips aiid work in a frame where the plus components of q and k are large. In the 
high-energy limit, the amplitude corresponding to Fig. is given in the rest frame of 
the proton by 



^27r5(A;o-go)T"(Aj 



-tX (-A;+A^(A)) (e*A 



(2) 



Here A = k — q, the field A is the Fourier transform of the external color field A, and e is 
the polarization vector of the produced gluon. The evaluation of the r.h. side of Eq. (H) in 
the high-energy limit shows that T"^ does not depend on A+ and A_ due to the softness 
of the external field. It is convenient to consider the impact parameter space amplitude 



iXqo 
2Ca 



dx^ tr 



(3) 



(with 9x = d/dx±), which is related to the amplitude in Eq. (0) by a Fourier transforma- 
tion in transverse space. The x_ dependence of A-^ is irrelevant in the high-energy limit. 
Here A-^ = A^T^ and are the generators of SU{Nc) in the adjoint representation, 
which we use throughout this paper; Ca = Nc- 

Resumming the gluon exchange to all orders means that the fast gluon created at 
the initial X99 vertex acquires a non-Abelian eikonal factor while travelling through the 
rest of the external field. This is illustrated in Fig. 0. 



X 
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Figure ^ The process x — (7 in an external color field. The dot symbolizes the initial 
Xgg vertex. 



The transition to the resummed amplitude of Fig. 1^ is realized by the substitution 



in Eq. (|), where 



Ufoo,:c+)i.Xl_) = Pexp 



-— / dx+A_{x+,x±] 



(5) 



The operator P denotes path ordering along x+. 
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From this, the total cross section for the scattering of the virtual photon off the color 
field target, i.e., the semiclassical (sc) leading-order result, can be derived. Using the 
identity 

dx+ U(to^,^^{x^) {d^A-^{x+, xx)) f/(^+,_oo)(a;±) = ( ) d^U(t,,-oc){x±) , (6) 
one finds 



where the index (0) stands for 'leading order'. Here 

W^^iy^) = U^{x^)U^\x^ + - 1 (8) 

and 

U^{xA-) = Uf^,-oo){^A-)- (9) 

In Eq. (|^, the summation over color indices is implicit. In the following it is always 
assumed that an appropriate averaging over the color fields underlying the basic quantity 
W-^ is performed (cf. 0). 

y.......... 

Figure |^: The partonic process X9 ~^ 9 leading order. 

The leading-order semiclassical gluon distribution is obtained by comparing this re- 
sult with a conventional partonic calculation, where the target is described by a gluon 
distribution. The relevant diagram is shown in Fig. 0, and the corresponding parton 
model (pm) cross section reads 

afl{x,Q') = ^xg^'\x,Q'). (10) 



Identifying the cross sections of Eqs. (0) and (p!OD, one obtains 

x9^'\x,Q') = / d\^\d^W^^^{0^)\\ (11) 

This result has been derived in identifying the scaling violations of F2 with the gluon 
distribution. Note the color factor in Eq- ( P^ due to the adjoint representation, 
which we use throughout the paper. As expected, this leading-order gluon distribution 
xg^^\x, Q"^) is constant for s ^ and shows no scaling violations beyond those induced 
by the explicit a^-factor on the r.h. side of Eq. ([IT]). In the value of lim^.^o ^^^'(a;, Q^) 
has been given in the stochastic vacuum model. 
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3 Semiclassical calculation at next-to-leading order 



The leading InQ^ calculation of the process x ^ 99 the semiclassical approach re- 
produces the conventional gluon-gluon splitting function, as shown in 0. Here, we need 
the complete next-to-leading order total cross section for the scattering of a x particle 
off an external color field. This calculation can be simplified significantly if one starts 
with a discussion of all diagrams contributing to the forward amplitude. The total cross 
section follows from the imaginary part of this amplitude. Furthermore, it is convenient 
to begin by considering the two-gluon exchange approximation. 

Clearly, the leading order diagram for the forward amplitude is simply the square of 
Fig. |1]. At next-to-leading order, all the diagrams in Fig. ^ have to be considered. 




Figure |^: Two-gluon exchange diagrams for the forward scattering of the scalar 'photon' 
X off a color field. The mirrored analogues of c), d) and e) and the analogue of graph a) 
with the two t channel gluons attached to different horizontal lines are not shown. 

In the high-energy limit, the diagrams d), e) (with their mirrored analogues) and 
diagram f) do not contribute. This can be understood intuitively by recalling that the 
gluon field is localized in a given region of space. Therefore, in the limit of infinite plus 
momentum, the right-moving gluonic degrees of freedom have no time for a virtual fluc- 
tuation between their first and second interaction with the external field (see Appendix ^ 
for a more technical argument). A related discussion in the case of particle radiation in 
high-energy scattering of external fields can be found in |TB| . 

What remains to be calculated is the imaginary part of diagram a) (and the cor- 
responding diagram where the two t channel gluons are attached to different lines), 
diagram b) and diagram c) (with its mirrored analogue). This corresponds to calculat- 
ing the cross section for x ^ 99 from the amplitude defined by the three diagrams in 
Fig. ^. Note that this is a significant simplification since the amplitude in Fig. ^ has no 
contributions arising from cutting diagrams d), e) and f) of Fig. ^ 

The next step is to resum the interaction with the external field to all orders, i.e., to 
repeat the step leading from Fig. |l| to Fig. ||. The result is illustrated in Fig. ||. Here, the 
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Figure |5|: The simplified amplitude for x ^ 99i which can be used for the calculation 
of the total cross section. 



first diagram corresponds to the fluctuation of the incoming 'photon' into a gg pair before 
the target and the subsequent scattering of the gluons off the color field, treated in the 
eikonal approximation. In an expansion in powers of the external field, the leading term 
reproduces the first two diagrams of Fig. ^ The second diagram of Fig. ^ corresponds to 
the creation of the two gluons in the space-time region of the external field, via a X999 
vertex. The two fast gluons then acquire non-Abelian eikonal phases while travelling 
through the rest of the color field. An expansion in powers of the external field generates 
the third diagram of Fig. |^. 




Figure ^: The relevant contributions to the amplitude in the next-to-leading order 
semiclassical calculation. The dot in the second diagram symbolizes the initial X999 
vertex. 



The amplitude of Fig. ^ can be calculated using the methods of [|] (see also 
Some details relevant to this particular process are discussed in Appendix |B[ The result 
can be given in the form 



(12) 

where a and 1 — a are the longitudinal momentum fractions carried by the two gluons, 
N"^ = a{l — a)Q'^ , k± and k'j_ are the transverse momenta of one of the two gluons before 
and after the interaction with the external field respectively, and W-^ is the Fourier 
transform of the function defined in Eq. (§). Summation over the indices i,j and over 
the color indices of W"^ is implicit. 

The full next-to-leading order semiclassical result is given by the sum of Eqs. (|^ 
and (0): 

aUx,Q')=ai'J{x,Q') + ail\x,Q'). (13) 



4 Parton model result at next-to-leading order 

Working in c? = 4 + e dimensions, the amplitude corresponding to Fig. |^, which includes 
a factor 1/(2 + e) for initial state gluon polarization, a factor 1/2 for identical final state 
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particles, and the color factor Ca = N^, reads 

Here Ad = A/i~^/^ and gd = g^i-'"^"^ are the rf-dimensional couplings, and s, t, u are the 
usual Mandelstam variables of a 2-to-2 scattering process. 





+ 




Figure |^ Conventional partonic amplitude for the process X9 ~^ 99- The crossed con- 
tribution belonging to the last of the three graphs and the ghost diagrams that have to 
be added for covariant polarization summation are not shown. 

From this squared amplitude, the total partonic cross section a^^\z,Q'^) for the 
process X9 ~^ 99 is obtained by standard methods. The variable z is defined hj z = 
Q'^/iQ'^ + Alternatively, it is given hj z = x/y, where y is the fraction of the target 
momentum carried by the struck gluon. Combining this with the leading order result, 
one can write 



a. 



pm 



1 dz 



X Z 



{ao{d}S{l-z) + a'^'\z,Q')) yg,{y) 



(15) 



where o"o{d} = 7rAy4 (cf. Eq. (|TU|)). Note that in the following we only calculate terms 
enhanced by powers of ln(l/x). Therefore virtual corrections, which affect only the end- 
point z = 1 and do not produce such terms, do not contribute. 

Working in the MS scheme, we renormalize the bare gluon distribution gb{x) accord- 
ing to 

9b{x) = g{x,fi'^) - ^ —Pgg{z) I- + -fE -InATvl giy^fi"^) . (16) 

ZTT Jx Z L C J 

The result can finally be given in the form 



a. 



pm 



C^0{d=4} 



X 



— \5{l-z) + — 

X Z \ ZTT 



^?(l/,/i'), (17) 



where the integrand is only complete in the region 2; < 1. In Eq. (|T^ Pggi^) (with z < 1) 
is the usual gluon-gluon splitting function. The coefficient function C^^{z), characteristic 
of the process under consideration, has been derived to be 



C. 



MS/ 



Pgg{z) In 



1- z 



HCa 
Qz{l~z) 
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5 Extracting the gluon distribution 



As in the leading-order case discussed in Sect. ||, the gluon distribution is extracted 
from the next-to-leading order semiclassical calculation by identifying the next-to-leading 
order total cross sections of the semiclassical and the parton model approach given by 
Eqs. and ([l^): asc{x,Q^) = cXp^x^Q^). 

We write the gluon distribution xg{x,fi'^) entering Eq. (IT 

.2 



as 



(19) 



where xg^'^\x, fi'^) is given by Eq. ( [Tl| ) and xg^^\x, fi'^) is a higher-order correction. In- 
serting Eq. (|19D into Eq. (I^), both and g^^^ are kept in the 5(1 — z) term, but only 
the leading-order distribution g^'^^ is kept in the as contribution. Now, identifying asc 
and CTpm, the leading order parts of both cross sections cancel and one finds 



xg^^\x, n"^] 



N'^Sii + 2kikj ^~ 



4(27r) 

— — / dz 

271 



^ m + kl 



(20) 



When evaluating the r.h. side of this equation, all terms that are not enhanced by ln(l/a;) 
will be dropped. In particular, we can use 



2C, 



and 



cT(A - — 



z 12 



(21) 



Note that in Eq. (|20|) we are only interested in the leading twist contribution of the 
semiclassical cross section which will be determined in Eq. (p2|). 

The integral involving the function W-^ is conveniently rewritten using the integra- 
tion variable z = Q'^/{Q'^ + JVP), where = s = k'l/{a{l — a)) is the invariant mass 
of the two gluons in the final state. We introduce a parameter fi' such that, for a soft 
hadronic scale A governing the behaviour of W"^, one has A^<c/i'^ <Q^. This allows us to 
decompose the integral according to (see Appendix y for details) 



4(27r) 



^-r^zr [da [ .^^^ , [d'^x_L [ d^kx_- 
In)' J J ail — a) J J 



N^6ii + 2kikj ^~ 



W^ik'^-k^] 



where the non-perturbative W'^ dependence is encoded in the function /: 



(22) 



f{k't) 



(27r)2|/2 7 (27r)24 



(23) 
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Here and in the following we disregard all terms suppressed by powers of A^//i'^. The 
motivation for writing the integral in the form given in Eq. (^) is the explicit separation 
of the InQ^ term, which is multiplied by the short- distance specific function d±W^^{0±). 

Note that we have not specified the lower bound of the z integrations in Eq. (^21). 
Clearly, the kinematical limit is z = x since the invariant mass of the produced gg 
pair can not be larger then the total center-of-mass energy available. However, no such 
bound appears explicitly in the semiclassical treatment since the classical color field 
behaves like an infinitely heavy target. The physical cutoff is provided by the breakdown 
of the semiclassical approximation for z ^ x, i.e., y ~ 1. 

Therefore, in order to obtain the leading logarithm, the second term on the r.h. side 
of Eq. (I22D has to be treated by applying the substitution 

^ dz ^ 1 _ 

— — ^ In-. (24) 

z X 

For the first term on the r.h. side of Eq. (0), only the short distance structure of the ex- 
ternal field matters. According to Eq. (p!TD , this short distance behaviour is characterized 
by the leading order gluon distribution yg^'^\y, fx^). Thus, we substitute 

2^/t^"|^/^'"^I^^^-(°^)I' ^ £^ In |^2.r^a.(/.^) ,,(")(,, /.^) (25) 

and assume that the behaviour of the phenomenological gluon distribution at z ^ x, 
i.e., y = x/ z ^ 1, correctly accounts for the region where the semiclassical treatment is 



no longer valid. In fact, when inserting Eq. ( p2D with the substitutions Eq. (p^ ) and (25) 
into Eq. (pO]) , the details of the large-y behaviour of yg^^\y,^'^) do not matter since 
the ln^(l/x) enhanced contributions from the semiclassical and the partonic calculations 
cancel. Thus, we obtain 



xg^^\x, /i^) = In — 

X 



- C^A ( In ^ + ^ ) (X, /.^) + - r dk'lfik"' 



(26) 



where the /i'^ dependence cancels between the two terms. Therefore, we can set /i'^ 
/i^ exp[ll/12] and write 



xg^\x^^x') = -(\n^-) / dk'lf{k'l). (27) 

TT V Xj Jo 

Note that the correction xg^^^x, /i'^) shows scaling violations consistent with the 
Altarelli-Parisi evolution at small x and a logarithmic small-x enhancement that is sen- 
sitive to the non-perturbative, large-size structure of the target. 

It would be interesting to evaluate Eq. (27) in the framework of the model of the 
stochastic vacuum of or of the large hadron model employed in ||^. Note in particular 
that, following 0, the large- Ai'c expression 

/ d^xi_ tr {w^^{y^)W^I{z±)) = QN^ [l - e^'^^/i _ g— 1 + g-afex-^x)^] (28) 
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can be used in Eq. (p3|). Here is the geometrical cross section of the target, the impact 
parameter dependence of the target thickness (which would be reflected in an impact pa- 
rameter dependence of the parameter a) is neglected, and averaging over all relevant color 
field configurations is assumed. With this model, the unintegrated, i.e., k'^ dependent, 
version of xg^^\x, fx'^) reproduces the recent result of Mueller (cf. Eq. (44) of |1^). The 
result presented in Eq. (|27|) extends the discussion given in by carefully matching 
the semiclassical and the parton model calculations. This gives rise to a precise definition 
of the cutoff of the k'j_ integration in Eq. (pTf ) and to the interpretation of the result as 
a correction to the leading order gluon distribution of Eq. (|r^). 

Equations (p3D and (pT]) can be significantly simplified (cf. Appendix for details) 
leading to the main result of our paper: 



It is easy to see how such a formula comes about: if the cutoff /i^ in Eq. (^) could be taken 
to infinity, the integration dk'l = d'^k'^/ir would give rise to the 5-function 6'^{y± — z±) 
(cf. Eq. (|23|)), and either the y± or the z± integration could be trivially performed. The 
remaining integration, say the y± integration, is now divergent at small y±, showing that 
the cutoff in Eq. (^) can, in fact, not be removed. The integrand of this divergent y± 
integration is precisely the one of Eq. ( p9D and, as can be seen by closer inspection of 
the relevant integrals (cf. Appendix 0), it is possible to translate the upper cutoff of the 
k'j_ integration into a lower cutoff of the y± integration. This is the origin of the compact 
formula in Eq. 



Given a specific model for the gluon fields of the target that allows for the calculation 
of the fundamental quantity W^^{y±), Eq. (|29|) can be used to improve the leading order 
semiclassical result of Eq. (0) . 

With the large- iVc expression of 0, 



-2 



/ dVtriy^^Jy^) = J dVtr {w,^^{y^)W.^I{y^)) = 2^1 [l - e'^^^i] , (30) 



one can see the qualitative agreement of Eq. (^9]) with Eqs. (49)-(51) of |T^. The progress 



of the present investigation is the careful matching of the semiclassical and the parton 
model treatment, giving rise to the precise definition of the cutoff of the y±_ integration 
in terms of the scale of the gluon distribution. 



6 Summary and conclusions 

The semiclassical approach has been successfully applied to different kinds of high-energy 
reactions. In DIS, the interaction of the energetic partons of the photon with the tar- 
get can be calculated in the eikonal approximation and is essentially given by a (non- 
perturbative) Wegner- Wilson loop, which measures an integral of the field strength of 
the target. 
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This approach, though hmited to low scales /i^, can predict input distributions for 
the DGLAP equation. Of course, in order to make numerical estimates of the parton 
densities, the evaluation of the Wegner- Wilson loop in a specific non-perturbative model 
is necessary. 

In this respect the gluon distribution is of particular interest as it dominates DIS 
at low X. To leading order, the distribution xg^^\x^iJ?) is a constant characterizing the 
averaged local field strength in the target. Therefore, a calculation at next-to-leading 
order is mandatory to obtain a nontrivial energy dependence of the gluon density. 

Using a scalar 'photon' (denoted by x), which couples directly to the gluon in a gauge 
invariant way, we derive the gluon density by matching the semiclassical and the parton 
model approach. At leading order, we have to equate the cross section for the transition 
X — ^ fif in a soft external field with the cross section of the process X9 ^ 9 given in the 
parton model. At next-to-leading order, the cross section for the transition x ^ (7(7 in a 
soft external field has to be equated with the parton model cross section of the process 
X9 ~^ 99- The ag correction of the gluon distribution shows a logarithmic enhancement, 
i.e., X g^^\x^ jj?) oc ln(l/x). 

In our approach we can not go beyond the leading-ln(l/x) approximation. This is 
a fundamental limitation of the semiclassical treatment and many related approaches, 
which has its origin in the artificial separation of the QCD dynamics into the soft degrees 
of freedom of the target and the high-energy modes of the projectile. Integrating over 
the fluctuations of the projectile, one has to drop the soft modes since their interaction 
can not be treated in the eikonal approximation. 

Parton distributions at next-to-leading order are scheme dependent. Our final result 
in Eq. ( ^9]) (together with the leading order contribution of Eq. ([TI|) ) provides the gluon 
density in the MS scheme. The short- distance cutoff in Eq. ( p9| ) is quantitatively related 
to the scale of the gluon distribution. This enables us to obtain numerical predictions 
for the gluon distribution at next-to-leading order in any non-perturbative approach 
describing the soft color field of the proton. 
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A Suppression of intermediate virtual fluctuations 
in the high-energy limit 

To illustrate the vanishing of diagrams d), e) and f) of Fig. ^ in the high-energy limit, 
consider first the simpler case where the horizontal gluon lines are replaced by scalar 
lines. The scalar version of diagram Fig. ^) with the naming of momenta used in this 
appendix is shown in Fig. ^. 

In the high-energy limit, the soft external field can not change the momentum plus 
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Figure |^: Scalar version of diagram Fig. |d). 

component of the fast right-moving particles essentially, ~ g+. Furthermore, the 
dependence of the external field vertex on the minus component of the transferred mo- 
mentum can be neglected. Thus, the minus component integrations in the diagram Fig. | 
take the form 

"^^-Z"^^- [k+k_ -kl + te] [/+/_ -11 + te] [{k-l)+{k-l)_ - {k-l)l + It] ■ ^^^^ 

Note first that this expression vanishes unless < /+ < A;+, in which case the two poles 
in the complex /_ plane lie on different sides of the integration contour. However, this 
means that both kj^ and {k — /)+ have to be positive, so that now both poles in the 
complex k^ plane are on the same side of the integration contour. This demonstrates 
that the integral / vanishes in any case. 

The same argument applies to the scalar analogues of diagrams Fig. ^) and f), since 
there the momentum k also flows through more than one propagator between the two 
external field vertices. This results in the presence of several poles in A;_, all on the same 
side of the integration contour, and therefore in the vanishing of the integral. 

All that has been said above applies to the gluonic case as well, as long as there is 
no additional minus component dependence introduced by the numerator factors. Such a 
dependence would prevent us from closing the integration contour. Working in Feynman 
gauge and applying the decomposition of the gluon propagator given in Appendix ^ to 
the lines with momentum /c, one can easily convince oneself that, in the high-energy 
limit, there is indeed no dependence on the minus components. 

Thus, we are justified in disregarding diagrams Fig. ^) - f) when calculating the 
forward scattering amplitude at next-to-leading order. 



B Details of the next-to-leading order semiclassical 
calculation 

The l.h. diagram of Fig. ^ contains three contributions: two of them describe the interac- 
tion of a single outgoing gluon with the color field, whereas in the third case both gluons 
interact with the field. Here we consider the latter case as an example, and outline some 
crucial steps of the evaluation of the corresponding amplitude which we call T^^. 

Working in Feynman gauge and denoting the momenta of the gluons before (after) 
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the interaction by k and p {k' and p'), respectively, the amphtude reads 



-^AW)^^^^P'^P)^H,„{p, k)^V^,{k', k)e'*^{k') , 

(32) 

where the indices i,j characterize the polarization of the gluons in the final state. The 
quantity H represents the X99 vertex, 

Hpa{p, k) = i\{{kp) gp„ - kpp„) , (33) 

while V{p' ,p) and V{k' , k), which will be explicitly given below, are the effective vertices 
for the scattering of the fast gluons off the external field. 

To simplify the calculation we exploit an alternative representation of the metric 
tensor appearing in the gluon propagators. For instance, the tensor g^P of the propagator 
—ig^^P Ip^ is decomposed according to 

with the light like vector n = {n+,n_,n^) = (0,2,0^). A possible choice of the po- 
larization vectors e(j)(p) in Eq. (p^), which in particular satisfy e(j)(p)p = 0, is given 
by 

e(,)(p)=fo,^^^,e«^), (35) 



P+ 

with the transverse basis e(i)_L = (1,0) and €(2)_l = (0, 1). 

In the high energy limit, the t channel exchange of gluons leads to an amplitude 
which is proportional to g+. As a consequence, in order to obtain the leading part of 
the amplitude, in Eq. (|3^) only the term containing the polarization vectors has to be 
considered. Because of gauge invariance of the X99 vertex the second term on the r.h. side 
of Eq. (0) vanishes. The contribution of the remaining two terms is at most proportional 
to {q+)^ and hence negligible. 

The effective vertices are now multiplied by the polarization vectors of two on-shell 
gluons; these vertices are governed by the non-Abelian eikonal factor defined in Eqs. (P) 
and d^). One finds 

e1^{p')V^p{p\p)e>l^^{p) = 2n5{p',-po)2po[U^ip'^-p^) - {27rY5'ip'^ - p^)]5,, , (36) 

and an analogous expression for the second vertex. The contribution proportional to 
S'^{p'^ — p_i_) subtracts the term in U-^ which contains no interaction, while 6ij indicates 
conservation of the gluon helicity. In Ref. the interested reader can find more details 
of the derivation of Eq. (^). 

Writing the loop integration in terms of light cone variables, d'^k = {l/2)dk+dk-d'^k±, 
and using the approximation 6{kQ — ko) ~ 25(A;'^ — /;;+) we arrive at 

%f = I - + P')h + '^hkj] (37) 

U\A-P±) - {27:f 6\p'^-p^)] [U^\k'^ - k^) - {2'nf6\k'^ - k, 
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The diagram on the r.h. side of Fig. |^ cancels the k"^ and p"^ term in the expression 
(Q^ + k'^ + p^) of Eq. (|37D . Now, the /c_ integration can be performed by closing the 
integration contour in the lower half of the complex fc_ plane. 

The resulting expression contains terms proportional to U-^U-^\ U'^, and U'^\ as 
well as a constant term. The two additional contributions, where only a single gluon 
interacts with the external field and which belong to the l.h. side of Fig. |^, contain terms 
proportional to f/"^ and f/"^^ and a constant term. If one adds these contributions, the 
color field dependence of the total amplitude Tij turns out to be 

U^{pl~P±)U''Hk'± - k^) - i27T)^6\p'^-p^)6\k'^ - k^) (38) 
= j rf2a;^e-*"^('=^+P'^)iy,^(A;l - ki_) . 

This provides us with the final result of the total amplitude, 

'^- = ^J '^'- ''^nCT' S ^'^-^'^^^^'^"^'^^^^^J^l - k.) . (39) 

The cross section in Eq. (|I^) is obtained by means of the standard formula for 
scattering off an external field, 

where the two identical particles in the final state require the factor 1/2 in front of the 
integral. 



C Evaluation of the integral with the function W'^ 

In this appendix we present some details of the derivation of Eq. (E^I). To introduce 
the integration variable z we exploit the relation k'l/{a{l - a)) = = Q^{1 - z)/z. 
Therefore the k'j_ integration appearing on the l.h. side of Eq. ( P^ can be replaced 
according to 

dk'f 2 dz 



Q' -V^ (41) 



a{l — a) ^ J z'^ 

where just the leading contribution at small z is kept. Concerning the lower bound of 
the z integration we refer to the discussion given in Sect. |[ 

Subsequently, it is convenient to divide the a integration of Eq. (^) in two parts 
by introducing an arbitrary parameter fi' which fulfills the condition <^ fJ'''^^Q^- 
To be specific we separate the symmetric and asymmetric gluon configurations in the 
integration by means of 



da= da+ 2 da, (42) 

J JtM'yNP Jo 
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where the symmetric (asymmetric) configurations give rise to the first (second) term on 
the r.h. side of Eq. (|22|) . For a given value of the symmetric contribution contains 
large transverse momenta k'j_ while the asymmetric term is entirely soft. 

To extract the leading twist of the symmetric hard part we use 



with C, 



[ m + k'f 



A 

X 



—kj_) + higher twist , (43) 
k'aiN%, 



The simple relations 



d'k±kaW^, (-k 



-' tj ,a'~^ ij ,b 
2 



Az{l - z){l + z'^ 
{2^Y djvt (0, 



(A^2 + ^ 



bah, and 

2 



a\2 



(44) 
(45) 



are important for the further evaluation of the hard contribution. In writing Eq. (|44|) we 
have anticipated the symmetric k'^ integration and employed the substitution k'Ji'^ 
k'lSab/'^- Finally, performing the a integration and keeping only the contribution that is 
dominant at small z, the first term on the r.h. side of Eq. ( ^2]) is obtained. 

In the case of the asymmetric soft contribution, which is independent of and 
therefore leading twist, the a integration is replaced by an integral over k'^. Since a is 
small we have Q'^a ~ k'lz, leading to the replacement 



da 



z 



dk 
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(46) 



The remaining step to obtain the second term on the r.h. side in Eq. (|22|) is the kj^ 
integration. Keeping in mind that for the soft term we can neglect N"^ this integration 
can easily be done with the aid of 



/ 



dkl^e^'^^y^ 

ki 



(47) 



D Simplification of the final expression 



which represents a particularly 



In this appendix, the details of the derivation of Eq. ([ 
simple form of our final result, are described. 

It is convenient to introduce a real parameter e > and to write the gluon distribution 
given by Eqs. (|23| ) and (^7|) in the form 



X(?«(x,/i2) = ^(ln-)Jo(0,/i'), 



where 



h{a,b) 



X 



d'k'Mk'l) , 



(48) 
(49) 
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(Py± f d'^z± 



-1 



and 



Furthermore, we have 



Jo(0, fi') = hm [J,(0, oo) - oo)] 



(50) 

(51) 

(52) 



The k'j_ integration in the definition of /^(O, oo) is easily performed giving a 5-function 
of y± — z±. The result is 



/.(0,oo) 



(27r)2yt 



Since 



/i(2/±, ^±) - C'?/^-2:_L for Ij/^I , < 1/A, 
where C is a constant, one can write 

^^^°'-)-2^27 + 4^/.<.,f ^^^-^-) + ^^^)• 



(53) 
(54) 

(55) 



The parameter r has been introduced to separate the small-distance from the large- 
distance part of the y± integration. 

In the k'j_ integration defining /^(/i', oo), the momentum variable k'j_ is always large. 
Therefore, the result is only sensitive to the small distance structure of h{y±,z±), and 
we can write 



J,(/i',oo)= / d'^k' 



d'^y± f d'^z± 



Cy±-z±e 



. y\A 



This integral can be calculated using standard methods: 



/,(/i', oo) 



C /i 



/-2e 



27r2 2e 



'r(i + e)r 



2 



rfi 



+ + 



(56) 
(57) 



When inserting Eqs. ( [55| ) and (0) in Eq. (|52D, the poles in e cancel. Setting the parameter 
r, which so far was arbitrary, to 



r(/i)^ 



,/2 



,2 ' 



the result of Eq. (p9D follows. Here 7e is Euler's constant. 



(5^ 
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